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Abstract. The crossing probabihty in the time direction, irt, is defined for an ofF- 
equilibrium reaction-diffusion system as the probabihty that the system of size L is 
still active at time t, in the finite-size scaling limit. Exact results are obtained for 
the diffusion-limited coalescence problem in 1 -I- 1 dimensions with periodic and free 
boundary conditions using empty interval methods, ttj is a scale-invariant universal 
function of an effective aspect ratio, L^/Dt, which is the natural scaling variable for 
this strongly anisotropic system. 



1. Introduction 

The study of crossing probabilities for standard percolation has been the subject of much 



interest during the last decade [1-16] (for a recent introductory review see reference [ll7| ). 

In two dimensions, the crossing probability may be defined as the probability vr 
to have at least one cluster joining two opposite edges of a rectangular-shaped finite 
system with length Ly and width Lj_. It turns out that, at the percolation threshold, 
in the finite-size scaling limit {L\\ oo, L± — oo, with r = L±/L\\ fixed), the crossing 
probability is a scale-invariant universal function^ of the aspect ratio r [l|, |^. 

Following the numerical work of Langlands et al [|^, Cardy was able to derive 
an exact expression for 7r(r). Using the relation between percolation and the g-state 



Potts model in the limit q ^ \ ||18| and boundary conformal field theory, he obtained 
the crossing probability between two non-overlapping segments on the edge of the half- 
plane at criticality. The corresponding result in the rectangular geometry was then 
obtained through a conformal mapping. The non-trivial scale invariance of 7r(r) is 
linked to the vanishing of the scaling dimension x{q) of a boundary condition changing 
operator of the Potts model in the percolation limit, q ^ 1. Let us notice that some of 
these results have been rigorously proven recently [^, [Tl| . A related problem concerns 
the number of incipient spanning clusters at criticality [12-16]. Exact formulas have 
been also obtained in this field through conformal and Coulomb-gas methods [|13|, O]. 



In a recent work [|19] the critical crossing probability was studied numerically for a 



strongly anisotropic system, namely, directed percolation in 1 -|- 1 dimensions. In this 
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case, the crossing probability in the time direction vr^ is also the probability that the 
system of size L remains active at time t. Anisotropic scaling |^0|, ^ then implies that 
the appropriate aspect ratio is r = L^/t, where z is the dynamical exponent. Here, too, 
it was found that, in the finite-size scaling limit, the critical crossing probability is a 
scale-invariant universal function of an effective aspect ratio which is the product of r 
by a non-universal constant. 

In the present work we continue the examination of critical crossing probabilities 
in strongly anisotropic systems by considering the case of diffusion-limited coalescence 
(DLC). This is one of the many actively studied off-equilibrium systems [22-26] which 
yields itself to an exact analysis [27-40]. We study the problem with periodic and 
free boundary conditions using the empty interval method, or interparticle distribution 
function method, which is reviewed in [^6], ^ |. 



The case of periodic boundary conditions is studied in section 2 using the standard 
empty interval method for a finite discrete system [33-35]. In section 3 we use a 
modification of the standard method to treat the problem with free boundary conditions. 
The results are discussed in section 3. 

2. Diffusion-limited coalescence with periodic boundary conditions 

We consider the time evolution of DLC on a one-dimensional lattice with L sites and 
periodic boundary conditions. Each site is in one of two states, either vacant or occupied 
by a particle A. The dynamics is governed by the following processes 

A0 

0A 

with the same rate D. When a particle jumps with rate D on a nearest-neighbour site 
which is already occupied, the two particles coalesce immediately on this site. Thus 
coagulation may occur either to the left or to the right. To simplify we assume that the 
L sites are occupied with probability one in the initial state at t = 0. As a consequence, 
the probability distribution of the particles A is translation invariant at later time t > 0. 



A0 0A (diffusion), A A (coagulation), (1) 



We study the time evolution of the system using the empty interval method |29 
Let the symbol • (o) denote an occupied (vacant) site. The probability for a given 
interval of length n to be empty at time t, 

n 

/„(t) = Prob(ST^^^), (2) 

is translation invariant on the periodic system with uniform initial conditions. Its time 
evolution involves the probability 

n n 

Fn{t) = Prob ( • b o • — o o" ) = Prob ( b o — ■ o o"*) (3) 

to have an empty interval of length n, either preceded or followed by an occupied site. 
Since 

n n+1 n 

Prob (b o ■ — o d'«) + Prob ( b o ■ — o o d' ) = Prob ( b o ■ — o d' ) , (4) 
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the following relation is obtained: 

F„(t) = J„(t)-4+i(t). (5) 
For n = 1, L — 1, the empty interval probability satisfies the master equation 

= 2D [Fn-i{t) - Fn{t)] = 2D [/„_i(t) - 2/„(t) + J„+i(t)] . (6) 

The gain terms correspond to processes in which a particle occupying the first site on 
the right (left) of an empty interval of length n — 1 jumps to the right (left) to diffuse 
or coalesce on the next site, thus leaving behind an empty interval of length n. The loss 
terms correspond to processes in which a nearby particle enters an empty interval of 
length n either from the left or from the right. The final form of the difference equation 
follows from (^. 

Equation (|^) has to be solved with the boundary conditions 

Jo(t) = l, /L(t) = 0. (7) 

The first one results from the expression of the site occupation probability, -Fo(^) = 
1 — /i (t) , the second follows from the fact that an initially non-empty system remains so 
at later time, since the coalescence process leaves at least one surviving particle. With 
a full lattice at t = 0, the initial condition corresponds to 

/n(0)=5„,0. (8) 

The master equation (|^) is solved through the Ansatz 

In{t) = J2Un)e-^'' (9) 

where (f)q{n) = UqSm{qn) +VqCOs{qn) when Uq = 8Dsin^(g/2) is non-vanishing. It takes 
the form (poi^) = an + b for the zero mode, c^o = 0, corresponding to the stationary 
state. The first boundary condition in (^, 

Io{t) = l = b + J2^qe-'^-'\ (10) 

leads to 6 = 1 and Vq = whereas the second, 

/^(t) = = l + aL + ^Uq sin(gL) e""^'* , (11) 

gives a = —1/L and sin(gL) = 0. Thus the empty interval probability can be written 

as 

4(t) = l-- + 5^CfcSin(^-^jexp^-8Dtsin2(^— jj . (12) 

The stationary state solution /„(oo) = 1—n/L corresponds to a single particle diffusing 
on the L sites so that, a site being occupied with probability 1/L, an interval of n sites 
is non-empty with probability with probability n/L. 
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According to (||), at t = we have: 

L-l 

^ Cfc sin 



k=l 



nkix 



L 



^-1, 
L 



Making use of the orthogonahty relation for the sines, 

L-l 



n=l 



sm 



/ nkir 



sm 



nln 



{k,l = l,L-l), 



in equation (|T3|), we obtain 

L 
2 



Cfc = — S2{k) 



where 



L 



Si{k) = n sin { — 

n=l ^ 
n=l ^ 



(13) 



(14) 



(15) 



-1^+^-cot — 
'2 \2L 



cot 



f k7i\ 
2l) 



so that finally 



exp 



^ ; , klT 

-SDtsin^ — 
2L 



(16) 



(17) 



The mean number of particles per site (or site occupation probability), p{t) = 1 — Ii{t) 
has the well-known long-time behaviour in the infinite system Here we 

are interested in the behaviour of the crossing probability in a system with aspect ratio 
r = L^/t with a dynamical exponent z = 2 for DLC The crossing probability in the 
time direction, Pt[L,t), is the probability that the system of size L is still active at 
time t. 

The probability that the system is in the stationary state, with the last particle on 
a given site, is equal to the probability lL_i{t) that the L — 1 other sites are empty. 
Since there are L possible choices for the occupied site, one obtains: 

L-l 



P,(L,t) = 1 -L/^_i(t) = 2X;(-l)'+'cos2(^) exp 

k=i ^ ^ 

In the finite-size scaling limit, this leads to the scale-invariant expression 

7r,(reff)= lim P,(L, t) = 2 V(-l)'=+i exp + 0(L- 



:i8) 



(19) 



where the crossing probability tt^ shown in figure [I| is a universal function of the effective 
aspect ratio Ves = r/D = L'^/Dt. 
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Figure 1. Variation of tlie scale- invariant crossing probability with the effective aspect 
ratio for periodic boundary conditions. The solid line corresponds to the asymptotic 
expression in equation 



3. Diffusion-limited coalescence with free boundary conditons 

3.1. Master equations for the empty interval probabilities 

With free boundary conditions a modified version of the empty interval method is needed 
to calculate the crossing probability. We define 



/ \ T-i 1 / -'^ n+l L 

Jm,n{t) = Prob (00---0 0---00 



(20) 



as the probability to have two disconnected empty intervals with sites 1 to m and n + l 
to L empty Its time evolution depends on the probabilities 



m+l n+l L 
O • O ■ ■ ■ O O 



Gm,n{t) = ( O O 

Hm,n{t) = Prob ( O O ■ ■ ■ O 

The different probabilities satisfy the relations 



m n L 
O ■ ■ ■ O O 



(21) 



Prob ( 



1 m n+l L , 

00---0 o---OOj= 

,1 m+l n.+l 

= Prob o o • ■ ■ o • o ■ 



L 

■ o o 



,1 m+l n+l 

Prob o o ■ ■ ■ o o o ■ 



L 

■ o o 



m n 
O WO 



L 
O O 



Prob ( o o 



m n 
O O O 



L 
O O 



(22) 



= Prob ( o o • 
so that we have 

Grm,n(t) Jm,n(t) ^m+l,n(^) j 

Hm,n(t) Jm,n(t) Jm,n—l(t) ■ (23) 

As before we assume that the L sites are occupied in the initial state. Thus the system 
contains at least one particle at later time. The condition of non-emptiness can be 
written as 

Jm,m(t) = 0, m = 0,L. (24) 
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One may notice that according to the definitions given in (pTD, Gm,m+i{t) = Hm,m+i{t) 
Jm,m+iit), in agreement with equations (ESf) and (El). 



When < m < n < L, the empty interval is indeed built of two disconnected parts 
and its probability satisfies the master equation 

^'^ ^ = D [Gm-l,n{t) + Hm,n+l{t) — Gm,n{t) — Hm,nit)] 

= D [Jm-l,nit) — 2Jm,n(t) + Jm+l,n(t) + J'm,n-l{t) ~ '^Jm,n{f) + Jm,n+lit)\ (25) 

where the gain terms correspond either to a particle at m jumping to the right or a 
particle at + 1 jumping to the left and the loss terms either to a particle at m + 1 
jumping to the left or a particle at n jumping to the right. When m = 0, there is a 
single empty interval from site + 1 to site L and the master equation reads 

^^^^ = D [i/o,n+l(t) - Ho,n{t)] = D [Jo,n-l(t) - 2Jo,n(t) + Jo,n+l(t)] • (26) 

In the same way, when n = L, we are left with a single empty interval from site 1 to 
site m and the corresponding probability evolves according to 

^^^^^ = D [Gm^iM - G„x(t)] = D [Jm-l,L{t) - '^JmAt) + J^+lA^)] ■ (27) 

Equations ( ^61) and (^Tf ) contain the same gain and loss terms as for the corresponding 
empty intervals in equation (pSj). They remain valid for n = L — 1 and m = 1, 
respectively, provided Jm,nif) satisfies the boundary condition 

^o,L(t) = 1 . (28) 

3.2. Solution of the eigenvalue problem 
Looking for the solutions under the form 

Jm,n (t) = ^ 0^ (m, n) e""* , (29) 

the master equations (p5| -p^) lead to an eigenvalue problem which has been discussed 



in details in reference [p4 |. 



Since < m < n < L, there is a total of L{L + l)/2 modes. According to (2i), 
(j)^^[m.,n) is an antisymmetric combination of eigenfunctions of the second difference 
operators involved in (p5|-p^). The problem is invariant under space reflection so that 
0(^(m, n) can be chosen as an eigenfunction of the space reflection operator 

P : (m, n) \—>- (L — n, L — m) . (30) 

Three types of solutions are obtained pl| : 

• The stationary solution 

JrriA'^) = 0o("^,n) = (31) 

which is an eigenstate of P with eigenvalue +1. Its expression follows from the 
fact that the zero mode eigenfunction is linear in m and n and has to satisfy 
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the boundary conditions ( p^ and (0). The first condition leads to the form 
(l)o{m,n) = c{n — m) and the second gives c = 1/L. It has also a simple physical 
interpretation: since a single particle remains in the stationary state, a site is 
occupied with probability 1/L. Thus the probability to have L — n + m empty 
sites, from 1 to m and from + 1 to L, is given by 1 — {L — n + m)/L. 
One may notice that with Jq^i{oo) = 1, the time-dependent part of Jm,nif) has to 
satisfy the boundary condition 



\L{t) - JoA^) = 

according to (p8|). 

The 2(L — 1) one-fermion excitations 



(32) 



1 


( mkn^ 




ri L , 


1 









sm 



mkTT 



2m 



sm 



nkiT 



(33) 



with k = 1, L — 1. These functions are antisymmetric eigenstates of P such that 
= ±(— 1)'^0^. They vanish when m = and n = L in agreement with (|32|). 
The excitations energies are given by 

. ( krc 



cjfc = 4:D sin 



(34) 



\2L^ 

Actually the odd eigenstate of P, (j)'^{m,n), is the combination of a one-fermion 
excitation with a zero mode. 

The (L — 1)(L — 2)/2 two-fermion excitations 



2 
L 



( mk7T\ , f nl-K 



sm 



sm 



— sm 



m/vr 



sm 



nkn 

~17 



(35) 



with \ <k <l < L — \. These antisymmetric two-particle states are eigenstates of 
P with eigenvalues (— l)'=+'+i and they satisfy the boundary condition (|32D . The 
corresponding eigenvalues read 



sm 



kTT\ 



— \ + sm — 



2L J 



ln\ 



- 



2LJ 



The solution satisfying the boundary conditions can be written as the expansion 



Jm,n (t) " 



n—m 



L 



L-l 
k=l 



E 



a2 d)t(m, n) 



L~2 L-l 

EE' 

k=i i=k+i 



ki (pki[m,n)e 



All the sites are occupied with probability one in the initial state, so that 

'-^m,n(0) ^m,0 ^n,L • 

Thus, for < m < n < L, we have: 



L-l 



L-2 L-l 



(36) 



(37) 



(3^ 



(39) 



k=l a=± 



k=l l=k+l 
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orthogonahty relations between the different eigenf unctions. Following reference |3| 
let us define surface and bulk scalar products for arbitrary functions / and g as 

L-l L-l 

{f\9)s = Yl /K ^) aim, L) + Y, /(O, n) g{0, n) , 



m=l 
L-l n-1 



n=l 



(/k)b = J] J]/Kn)(7(m,n). (40) 

n=2 m=l 

It turns out that the one-fermion eigenfunctions in (|33|) are orthogonal for the surface 
scalar product whereas the two-fermion eigenfunctions in (|35|) are orthogonal for the 
bulk one: 



l>)s =^kk'^af3 {k,k' = 1,L -1; a,f3 = ±) , 
(0fc/|0fc'r)b = hk' 5iv {0<k <l < L; <k' <l' < L) 
Thus, taking appropriate scalar products in (|39D, one obtains: 



(41) 



-(0^l0o)s 

I0o)s 



2S,{k) S2{k) 
S2{k) 



l-f-l)'^ 



l + fkn 

cot — 

2VL \2L 



L-l 

bkl = -(0H|0o)b - XlXl 
(-1)''- 



( kn 
cot I 



\2L^ 

(0fed0o)b 



k'=l a=± 



-IV k'n\ fln\ 

^ cot — cot — 

2L \2L) \2L) 



The relation 



L-l 



^0^,(m,n)(0^,|0o)s = <Po{m,n) 



(42) 



(43) 



k'=l 



has been used in the calculation of bki- The final form of Jm,n{t) satisfying the initial 
and boundary conditions follows from (|37D and (0). 



3.3. Crossing probability 

Since Jn-i,nif) gives the probability that all the sites are empty, except site n which is 
occupied by the last particle, the probability that the system is still active at time t is 
given by: 

L 

Pi(L,t) = l-^ J„_i,„(t) 



n=l 

L-l 



k=l 



5^a^X^0^(n-l,n) 



a=ib n=l 



L-2 L-l 

-EE". 

k=i i=k+i 



'kl 



y^^(t)ki{n-l,n) 



n=l 



. (44) 



Scale-invariant universal crossing probability 



9 



Straightforward but lengthy calculations lead to the following results for the different 
sums over n: 



5^0+(n-l,n)=O 

n=l 
L 

^0^(n- l,n) = - 

n=l 



2 ' cot — ^ 



n=l 

Thus we obtain 

. L-l 



>kl 



'n — 1, n) 



sin(^) sin(^) 



sm 



(k~l)-K 
2L 



sm 



{k+l)n 
2L 



k=l 
k odd 

L-2 L-l 



2L y 



/7r\ 



2L J 



fc=l l=k+l sm 2L 2L 



2L 



sm — - +sm — - 



TT 



2L 



(45) 



• (46) 



k+l odd 

In the finite-size scaling limit, the crossing probability displays the scale- invariant 
dependence on the effective aspect ratio rgfr = L'^/Dt and reads 

7rt(reff)= lim Pt{L,t) 

L,t—*oo 
r fixed 



k=l 
k odd 



exp 



k^TX^ 



32 ^ ^ f-ll 



fc=l l=k+l 

k+l odd 



exp 



(fc^ + /^)7r 

rcs 



+ 0(1-'). (47) 



The rapid convergence to the scaling limit is shown in figure ^ 




Figure 2. Variation of the scale-invariant crossing probability with the effective 
aspect ratio for free boundary conditions. The solid line corresponds to the asymptotic 
expression in equation (^). 
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4. Discussion 

We have studied DLC at the critical point where the particle density and other quantities 
display power laws in the thermodynamic limit. The problem can be made off-critical 
by introducing a birth process with rate A, corresponding to the baclc-reaction of the 
coagulation process in (0). The system being strongly anisotropic, under a change of 
the length scale by a factor b, the length transforms as L' = L/b whereas the time 
transformation, t' = t/b^, involves the anisotropy or dynamical exponent z For a 
scale-invariant crossing probability one obtains 

P,(L,t,A) = P,(^P^,6i/^A^ , (48) 

where u is the exponent of the correlation length, ^ = i A~^, which diverges at the 
critical point, A = 0. Taking b = A~'^ leads to 

t. A) = ( A., 1^ = / 1 ^ , (49) 

where we introduced the relaxation time, r = fA'^". Both the correlation length 
and the relaxation time contain a non-universal amplitude, ^ and f, respectively. In 
equation (^), f{x,y) is a universal scahng function of its dimensionless variables, 
X = L/C, and y = t/r. The finite-size scaling limit at criticality amounts to take 
A = and 6 = L in (ffSI), which gives 



Pt{L, t, 0) = Pt (^1, 0^ = TTi (^c . (50) 

Thus the crossing probability ttj is a scale-invariant universal function of the effective 
aspect ratio cr. The non-universal amplitude c depends on the choice of the length and 
time units. It can be expressed as a function of the non-universal correlation length and 
relaxation time amplitudes by comparing (^) to (^9|). Since r appears through the 



dimensionless ratio x^/y, one obtains 

cr^^^ir (51) 

t/T 

and = 7^/^^ which is equal to in our case. 

We have shown on the example of DLC that the crossing probability vr^ is a scale- 
invariant function of the effective aspect ratio for different types of boundary conditions. 
This function is expected it to be universal as in the case of directed percolation |T9|. An 



indication of the universality of tt^ can be found in where a birth process of the form 
A0A ^ AAA was added to (|l]). The problem stays in the same universality class 
and remains exactly solvable through the empty interval method for all values of A. It 
turns out that Init) is only modified by terms of higher order in 1/L which disappear in 
the finite-size scaling limit. One could also check the universality of tt^ on the diffusion- 
annihilation problem, which has been shown to belong to the same universality class as 
DLC, through a similarity transformation in the quantum Hamiltonian formulation of 
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the master equation [38-40]. Another possibihty would be to verify the independence 
on the initial conditions, provided a finite fraction of the sites are occupied in the initial 
state. 

Finally, let us mention that a recent generalization of local scale invariance for 
strongly anisotropic systems |^T], ^ leaves the hope for directly obtaining the crossing 
probability formulas in a given universality class, like in isotropic systems. 
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